

* FOR FIGURE 6

* SIMULATION RESULTS


*  COPY chg_energy_co2 parameter FROM results/simulationsresults_tomerge/merged.gdx 
save "simulation_results.dta", replace 
  
rename var1 spec
rename var2 region


  
* FIG6
twoway (lfit co2elastconsmrio logpci if spec =="NH_resprod_06_fitteddem_yes_ESTTHETA"  & region !="LUX" ,  lcolor(black) lpattern(solid) ) /// 
 (scatter co2elastconsmrio logpci if spec =="NH_resprod_06_fitteddem_yes_ESTTHETA"  & region !="LUX" ,  mcolor(black) msize(small) msymbol(x)) ///
 (lfit co2elastconsmrio logpci if spec =="NH_CLMquad_resprod_06_fitteddem_yes_ESTTHETA"  & region !="LUX" ,  lcolor(red) lpattern(solid) ) ///
 (scatter co2elastconsmrio logpci if spec =="NH_CLMquad_resprod_06_fitteddem_yes_ESTTHETA"  & region !="LUX"  & region !="QAT", mlabel(region) ///
   mlabsize(vsmall) mlabcolor(red) mcolor(red) msize(small) msymbol(+)) ///
 , legend(order(2 "CRIE" 1 "" 4 "NH CES quadratic" 3 "" )  cols(2) ring(0) position(2)  ) xlabel(#7) ylabel(#7) ///
 ytitle(Income elasticity of total CO2 content of consumption) xtitle(Log per capita income)  yline(1, lcolor(black))  graphregion(color(white))
graph export "..\figures\FIG6.pdf", as(pdf)  replace
 graph export "..\figures\elast_mrio_comparison_2ndrev_red.eps", as(eps) preview(on) replace

 
 
 
 
 
